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1 Introduction 

The electroencephalographic (EEG) or magneto- 
encephalographic (MEG) inverse problem is to es¬ 
timate the current dipole sources underlying meas¬ 
ured electric potentials or magnetic fields of the 
brain outside the head. The estimation of multiple 
dipole parameters is known as a difficult nonlinear 
optimization problem due to local minima in the 
cost function. In order to reduce the computational 
complexity, we tried to find independent active di¬ 
pole sources from sensor signals using independent 
component analysis (ICA) and to fit the each source 
to a single dipole. The process of the source separa¬ 
tion makes a multiple dipole problem into several 
single dipole problems. In this paper, we proposed 
ICA dipole fitting method and tested the method 
through computer simulations. 

2 Methods 
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If s^. ’s are statistically independent, each a qi is an 
independent brain topography due to each dipole q z . 
We can find 2i qi ’s by applying ICA to the measured 
EEG or MEG. 

2.2 Independent component analysis (ICA) 

ICA is an effective technique for finding statistically 
independent sources and independent basis vectors 
spanning the signal subspaces. ICA dipole fitting 
method assumes that the dipoles are fixed in un¬ 
known location and orientation and the strength of 
each dipole independently changes with time. We 
used infomax ICA algorithm [2]. Our ICA dipole 
fitting method is comprised of the following steps: 


2.1 Forward Problem 


We assume simultaneous recordings at m sensors 
for n time instances by p current dipole sources. 
The spatio-temporal data matrix is expressed as[l] 

F = [G(r ql ) ... G(r^)Iqf q T J (1) 


where mX 3 matrix G(r y/ ) is the gain matrix that 
maps a dipole at r qi into a set of measurements and 
3 X n matrix q f is the time course of the moment of 
a dipole at r qi . If we assume all dipole sources have 
fixed orientations during measurement, the data ma¬ 
trix can be expressed as 


F = [G(r„) G(r„) 


ki 

i- 

o 

-1 

C/3 

' • • ^ 

_1 

i 

o 


-1 

& 

C/3 

_1 


( 2 ) 


where q. = u .s .. u . is a unit norm dipole orien¬ 
tation vector and 1X n vector s . is a time course of 
the dipole strength of n time samples. 


Step 1. Decompose the measured data into signal 
and noise subspaces by an eigen-decomposition 
(ED). 

R r =UAU5+U/X (4) 

The number of dipole sources is assumed to be the 
same as the number of significant eigenvalues. 

Step 2. Project the measured data into the signal 
subspaces. 

F s = U^F (5) 

The process reduces the dimensionality of the data 
and improves the signal to noise ratio (SNR). 

Step 3. Apply ICA to the projected data. 

S = WF S (6) 

Find W and the independent time course of the di¬ 
pole strength, S by ICA. 

Step 4. Calculate independent topographies. 

A - U S W _1 (7) 

Each column vector of the matrix A represents an 
independent topographic map due to a dipole source. 



2.3 


Source localization 


1.5 


Strength of the three dipoles 


For each independent topographic map, we can es¬ 
timate parameters of a single dipole using optimiza¬ 
tion techniques. We used downhill simplex method 
for minimization of the following cost function. 


C(rJ = 


— G(r • )[g r (r • )G(r •)]"' G r (r • )a, 
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A dipole location is searched by the simplex method 
and the dipole orientation is calculated as the linear 
least square solution. By considering the topography 
of each dipole separately, we search for only one 
dipole and the search space becomes dramatically 
reduced. So we can search for the nonlinear least 
square solution without additional computationally 
complex global optimization techniques. 


2.4 Simulations 


We performed computer simulations to test the ef¬ 
fectiveness of the proposed method in noisy envi¬ 
ronment. MEG data were simulated for the case of 
three dipoles with fixed locations and orientations. 
All the dipoles were randomly placed with the same 
eccentricity. And the orientations of the three di¬ 
poles were fixed to be (1, 0, 0), (0, 0, 1) and (0, 0, 1) 
respectively. The time course of the each dipole 
strength was assumed as shown in Fig. 1. 120 mag¬ 
netometers were uniformly positioned over the sur¬ 
face of an imaginary spherical head with radius 
10cm. A total of 1000 time samples were generated 
and were corrupted by additive white Gaussian noise. 

3 Results 



-" 1.5 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 

0 100 200 300 400 500 600 700 800 900 1000 

sample 

Figure 1: Time series of strengths of three dipoles 
assumed in the computer simulations. 



Figure 2: Mean dipole localization error. 


Table 1 and Fig. 2 show the simulation results. The 
estimation errors for dipole locations were less than 
2mm. The localization error was larger when a di¬ 
pole was placed closer to the center of head sphere. 

4 Discussion 

The simulation results showed the ability of the pro¬ 
posed ICA dipole fitting method. Without global 
optimization techniques, the proposed method was 
able to solve the problem of multiple dipole estima¬ 
tion. The results of our proposed method were better 
than those of using simulated annealing [3]. But the 
proposed approach has some problems. ICA need 
enough samples for learning and is not robust to 
noise. 
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Table 1: Dipole localization error ( mean ±std [mm] ) 



Eccentricity 

0.50 

0.65 

0.80 

SNR 

OdB 

1.24 ±0.70 

0.64 ± 0.26 

0.30 ±0.16 

6dB 

0.75 ± 0.44 

0.58 ±0.30 

0.28 ±0.05 

12dB 

0.67 ±0.48 

0.54 ± 0.24 

0.23 ±0.05 

18dB 

0.63 ± 0.45 

0.52 ± 0.24 

0.20 ±0.03 

24dB 

0.65 ± 0.47 

0.53 ±0.22 

0.22 ± 0.02 
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